global root_dir = "`1'"

include "$root_dir/code/config/config.do"

cap noi log using ${log_dir}/table_A9_skill_premium.log, replace name(tabl)

capture noi {

	qui do ${code_dir}/config/tabletools.do

	use ${final_dir}/regression_dataset_from1970_tfacit1.dta, clear
	tab year, gen(YY_)
	gen clustering_at = "lse_id"

	estimates clear
	clonevar LSW = lswMPm_1995_a
	clonevar HSW = hswMPm_1995_a
	clonevar GDPGAP = lngdpgap_1995_a
	clonevar VAEMP = vaempMPm_1995_a
	clonevar GDPPC = gdppcMPm_1995_a
	clonevar stockown = k${depvar}_${ttt} 
	clonevar stockownzero = k${depvar}_${ttt}0  
	clonevar stockother = kNOT_${depvar}_${ttt} 
	clonevar stockotherzero = kNOT_${depvar}_${ttt}0
	clonevar spilloversown = spill${depvar}${ttt}_1995_a
	clonevar spilloversother = spillN${depvar}${ttt}_1995_a
	clonevar spilloversownzero = spill${depvar}${ttt}_1995_a0
	clonevar spilloversotherzero = spillN${depvar}${ttt}_1995_a0
	sort lse_id year
	bys lse_id : egen _total_${depvar}_${ttt}_1995 = sum(${depvar}_${ttt}) if year>=1995+2 & year <= 2009+2
	bys lse_id : egen total_${depvar}_${ttt}_1995 = max(_total_${depvar}_${ttt}_1995)
	drop _total_${depvar}_${ttt}_1995
	egen yearctry = group(year country_shr_1995) if year <= 2009
	egen yearindustry = group(year industry) if year <= 2009
	egen yearctryindustry = group(year industry country_shr_1995) if year <= 2009	

	*1) F + IY
	gen LSWHSW = LSW - HSW
	ppmlhdfe F2.${depvar}_${ttt} LSWHSW GDPGAP stockown stockownzero stockother stockotherzero spilloversown spilloversownzero spilloversother spilloversotherzero  if year>=1995 & maxweight_1995 < 1 & missing_weights_1995==0 & missing_spill_weights_1995 == 0 & total_${depvar}_${ttt}_1995>0, vce(cluster lse_id) absorb(lse_id yearindustry)
	hasvar spill* stock*, local(stockspill)
	estadd local f "\yes"
	estadd local iy "\yes"
	estadd local cy "\no"
	estadd local obs "{\num{`e(N)'}}"
	estadd local firms "{\num{`e(N_clust)'}}"
	estimates store col1995_1
	*2) VAEMP | F + IY
	ppmlhdfe F2.${depvar}_${ttt} LSWHSW GDPGAP VAEMP stockown stockownzero stockother stockotherzero spilloversown spilloversownzero spilloversother spilloversotherzero  if year>=1995 & maxweight_1995 < 1 & missing_weights_1995==0 & missing_spill_weights_1995 == 0 & total_${depvar}_${ttt}_1995>0, vce(cluster lse_id) absorb(lse_id yearindustry)
	hasvar spill* stock*, local(stockspill)
	estadd local f "\yes"
	estadd local iy "\yes"
	estadd local cy "\no"
	estadd local obs "{\num{`e(N)'}}"
	estadd local firms "{\num{`e(N_clust)'}}"
	estimates store col1995_2
	*3) GDPPC | F + IY
	ppmlhdfe F2.${depvar}_${ttt} LSWHSW GDPGAP GDPPC stockown stockownzero stockother stockotherzero spilloversown spilloversownzero spilloversother spilloversotherzero  if year>=1995 & maxweight_1995 < 1 & missing_weights_1995==0 & missing_spill_weights_1995 == 0 & total_${depvar}_${ttt}_1995>0, vce(cluster lse_id) absorb(lse_id yearindustry)
	hasvar spill* stock*, local(stockspill)
	estadd local f "\yes"
	estadd local iy "\yes"
	estadd local cy "\no"
	estadd local obs "{\num{`e(N)'}}"
	estadd local firms "{\num{`e(N_clust)'}}"
	estimates store col1995_3
	*4) F + IY + CY
	ppmlhdfe F2.${depvar}_${ttt} LSWHSW GDPGAP stockown stockownzero stockother stockotherzero spilloversown spilloversownzero spilloversother spilloversotherzero  if year>=1995 & maxweight_1995 < 1 & missing_weights_1995==0 & missing_spill_weights_1995 == 0 & total_${depvar}_${ttt}_1995>0, vce(cluster lse_id) absorb(lse_id yearindustry yearctry)
	hasvar spill* stock*, local(stockspill)
	estadd local f "\yes"
	estadd local iy "\yes"
	estadd local cy "\yes"
	estadd local obs "{\num{`e(N)'}}"
	estadd local firms "{\num{`e(N_clust)'}}"
	estimates store col1995_4
	*5) VAEMP | F + IY + CY
	ppmlhdfe F2.${depvar}_${ttt} LSWHSW GDPGAP VAEMP stockown stockownzero stockother stockotherzero spilloversown spilloversownzero spilloversother spilloversotherzero  if year>=1995 & maxweight_1995 < 1 & missing_weights_1995==0 & missing_spill_weights_1995 == 0 & total_${depvar}_${ttt}_1995>0, vce(cluster lse_id) absorb(lse_id yearindustry yearctry)
	hasvar spill* stock*, local(stockspill)
	estadd local f "\yes"
	estadd local iy "\yes"
	estadd local cy "\yes"
	estadd local obs "{\num{`e(N)'}}"
	estadd local firms "{\num{`e(N_clust)'}}"
	estimates store col1995_5
	*6) GDPPC | F + IY + CY
	ppmlhdfe F2.${depvar}_${ttt} LSWHSW GDPGAP GDPPC stockown stockownzero stockother stockotherzero spilloversown spilloversownzero spilloversother spilloversotherzero  if year>=1995 & maxweight_1995 < 1 & missing_weights_1995==0 & missing_spill_weights_1995 == 0 & total_${depvar}_${ttt}_1995>0, vce(cluster lse_id) absorb(lse_id yearindustry yearctry)
	hasvar spill* stock*, local(stockspill)
	estadd local f "\yes"
	estadd local iy "\yes"
	estadd local cy "\yes"
	estadd local obs "{\num{`e(N)'}}"
	estadd local firms "{\num{`e(N_clust)'}}"
	estimates store col1995_6
	drop LSWHSW LSW HSW GDPGAP VAEMP GDPPC

	clonevar LSW_foreign = lswMPm_shr4_foreign_1995_a
	clonevar HSW_foreign = hswMPm_shr4_foreign_1995_a
	gen LSWHSW_foreign = LSW_foreign - HSW_foreign	
	clonevar GDPgap_foreign = lngdpgap_shr_foreign_1995_a
	clonevar GDPPC_foreign = gdppcMPm_shr4_foreign_1995_a
	clonevar VAEMP_foreign = vaempMPm_shr4_foreign_1995_a

	*7) F + IY + CY | Foreign
	ppmlhdfe F2.${depvar}_${ttt} LSWHSW_foreign GDPgap_foreign stockown stockownzero stockother stockotherzero spilloversown spilloversownzero spilloversother spilloversotherzero if year>=1995 & maxweight_1995 < 1 & missing_weights_1995==0 & missing_spill_weights_1995 == 0 & total_${depvar}_${ttt}_1995>0, vce(cluster lse_id) absorb(lse_id yearindustry yearctry)
	hasvar spill* stock*, local(stockspill)
	estadd local f "\yes"
	estadd local iy "\yes"
	estadd local cy "\yes"
	estadd local obs "{\num{`e(N)'}}"
	estadd local firms "{\num{`e(N_clust)'}}"
	estimates store col1995_7
	*8) VAEMP | F + IY + CY | Foreign
	ppmlhdfe F2.${depvar}_${ttt} LSWHSW_foreign GDPgap_foreign VAEMP_foreign stockown stockownzero stockother stockotherzero spilloversown spilloversownzero spilloversother spilloversotherzero  if year>=1995 & maxweight_1995 < 1 & missing_weights_1995==0 & missing_spill_weights_1995 == 0 & total_${depvar}_${ttt}_1995>0, vce(cluster lse_id) absorb(lse_id yearindustry yearctry)
	hasvar spill* stock*, local(stockspill)
	estadd local f "\yes"
	estadd local iy "\yes"
	estadd local cy "\yes"
	estadd local obs "{\num{`e(N)'}}"
	estadd local firms "{\num{`e(N_clust)'}}"
	estimates store col1995_8
	*9) GDPPC | F + IY + CY | Foreign
	ppmlhdfe F2.${depvar}_${ttt} LSWHSW_foreign GDPgap_foreign GDPPC_foreign stockown stockownzero stockother stockotherzero spilloversown spilloversownzero spilloversother spilloversotherzero if year>=1995 & maxweight_1995 < 1 & missing_weights_1995==0 & missing_spill_weights_1995 == 0 & total_${depvar}_${ttt}_1995>0, vce(cluster lse_id) absorb(lse_id yearindustry yearctry)
	hasvar spill* stock*, local(stockspill)
	estadd local f "\yes"
	estadd local iy "\yes"
	estadd local cy "\yes"
	estadd local obs "{\num{`e(N)'}}"
	estadd local firms "{\num{`e(N_clust)'}}"
	estimates store col1995_9
	drop LSWHSW_foreign GDPgap_foreign VAEMP_foreign HSW_foreign LSW_foreign GDPPC_foreign

	drop stockown stockownzero stockother stockotherzero
	drop spilloversown spilloversownzero spilloversother spilloversotherzero


	setlabels

	dhoztab * using ${tab_dir}/appendix/Table_A9_skill_premium.tex, b(%10.2f) se(%10.2f) nogaps replace numbers nonotes nolines nomtitles sfmt(a1) scalars("has_stockspill \noalign{\medskip} Stocks and spillovers" "f Firm fixed effects" "iy Industry \stimes year fixed effects" "cy Country \stimes year fixed effects" "obs \noalign{\medskip} Observations"  "firms Number of firms") label mgroups("Domestic and foreign" "Foreign", pattern(1 0 0 0 0 0 1 0 0) span prefix(\multicolumn{@span}{c}{) suffix(}) erepeat(\cmidrule(lr){@span})) drop(*zero _cons spill* stock*) noobs rename(LSW_foreign LSW HSW_foreign HSW GDPgap_foreign GDPGAP VAEMP_foreign VAEMP GDPPC_foreign GDPPC LSWHSW_foreign LSWHSW) ///
		depvar("Auto95") ///
		notes("This table shows the effect of the skillpremium on automation innovations. All columns include firm and industry-year fixed effects. Columns 4--9 add country-year fixed effects. Columns 7--9 compute the normalized foreign (log) inverse skill premium as the difference between the normalized (log) foreign low-skill wages and the normalized (log) foreign high-skill wages as defined in the text. In these columns, GDP gap, GDP per capita and labor productivity also correspond to their normalized foreign values. Standard errors are clustered at the firm-level and reported in parentheses.")
	*again, brackets
}
if _rc == 0 {
    display "Execution finished successfully."
}
else {
    display "Execution finished with errors."
}